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PREDICTION OF STRESSES IN AIRCRAFT PANELS 
SUBJECTED TO ACOUSTIC FORCES 

By 

Chuh Mei* 

SUMMARY 

This report summarized the progress and accomplishments performed under 
NASA/Langley Research Center Master Agreement NAS1-17993, Task Assignment 
No. 22, entitled "Prediction of Stresses in Aircraft Panels Subjected to 
Acoustic Forces," for the period October 1, 1985 to December 15, 1987. The 
primary effort of this task is the development of analytical methods for 
prediction of stresses in aircraft panels subjected to acoustic forces. The 
progress and accomplishments of various activities are discussed first. 

Then, publications, presentations and thesis are followed. 


♦Professor, Department of Mechanical Engineering and Mechanics, Old Dominion 
University, Norfolk, Virginia 23529-0247. 



PROGRESS AND ACCOMPLISHMENTS 

Seven activities were supported under the Task Assignment No. 22. The 
progress and accomplishments of these activities are sunmarized in the 
following. 

ACTIVITY 1: Analysis of Panel Random Response Using Laminate Theories of 

Various Degree of Detail. 

Two laminated plate theories are used: the classical laminated plate 

theory (CPT) and the first-order shear deformation plate theory (FSDT by 
Yang-Norri s-St avsky and Reissner-Mindl in) . Large deflection effects are 
also considered in the formulation. Results showed that the influence of 
transverse shear deformation on random response is considerable for moder- 
ately thick laminates (a/h < 20), and the small deflection linear theory 
with shear deformation included would give accurate predictions. For thin 
1 aninates (a/h > 50), the linear and nonlinear results agree well at low 
pressure levels, but disagree at high levels. Thus, the large deflection 
theory with shear deformation neglected should be used for thin plates at 
high levels of excitation. Therefore, for a particular value of a/h, one of 
the three simpler theories can be chosen that will provide accuracy equal to 
the more cimbersome theory which includes both shear and large deflection 
effects (see Publication P-1 and Report R-l). 

ACTIVITY 2: Effects of Nonlinear Damping on Large Amplitude Response of 

Beams and Plates. 

The teststrain spectral data of aircraft panels have repeatedly shown 
the broadening behavior of the peaks when the panels were excited at high 
pressure levels. It has demonstrated analytically, the first time, that 
nonlinear damping terms (c^q q, and C 2 qq ) will give the broadening 
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behavior of the strain peaks. Analytical results also showed that it would 
predict more realistic frequency of vibrations and RMS deflection with non- 
linear damping considered in the analysis. 

Further analytical efforts are needed to better understand the effects 
of nonlinear damping on panel response. Experiments are also needed to 
determine the nonlinear damping coefficients ( c ^ , c 2 etc.) for beam and 
plate structures . Other possible nonlinear damping terms should also be 
considered (Publications P-2 to P-5* and Report R-2). 

ACTIVITY 3: Development of Finite Element Large Deflection Random Analysis 

with a Beam Element and a Rectangular Plate Element. 

A literature search on this subject was conducted in June 1986. The 
development of finite element formulation using single-mode and multiple- 
mode for large deflection random response analysis of beams and plates has 
been completed. Finite element single-mode results agreed well with solu- 
tions using the Fokker-Pl anck-Kolmogorov (exact solution) equation and the 
equivalent linearization (EL) methods. Finite element (FE) multiple-mode 
beam results are shown in Figs. 1 and 2. RMS maximum strains and equivalent 
linear frequencies vs. sound spectrum level for simply supported and clamped 
beams using 3 modes in the analyses are shown in Figs. 1 and 2, respec- 
tively. Th e EL 3-mode solutions are also given for comparison. The finite 
element results by neglecting the axial displacement and the axial inertia 
degrees-of-freedom are denoted by FE* in Figs. 1 and 2. The FE* results 
agree well with the EL solutions, since both approaches are formulated in 


♦Publication P-4 received the best paper award presented at the AIAA 10th 
Aeronautics Conference. 
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RMS MICRO-STRAIN 



RE 2 x 10 -5 Pa 

Fig. 1. RMS maximum strain vs. sound spectrum level for 

simply supported and clamped beams (3-mode and c = 0.01). 
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Fig. 2. Equivalent linear frequency vs. sound spectrum level 

for simply supported and clamped beams (3-mode and ? =0.01). 
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terms of the transverse displacement only. 

Work is continued on the development of finite element large deflection 
random response of built-up structures (Publications P-6 and P-7). 

ACTIVITY 5: Comparison of Time Domain and Equivalent Linearization 

Sol utions. 

A comparison of maximum deflection spectra using the time domain and 
the EL methods for a 12 x 12 x 0.040 in. clamped, (0, ± 45, 90) $ composite 
plate is shown in Fig. 3. The time domain solution in Fig. 3 is obtained by 
Dr. E. Thomas Moyer of George Washington University. At 70 dB sound spec- 
trum level, the two methods predict very similar responses. At 110 dB, the 
EL method predicts essentially the same equivalent linear frequency as the 
time domain simulation; however, the time domain approach predicts the band- 
width of the response peak about 2 to 3 times broader than that of the EL 
method. The broader spectral response yields a 10% higher value of the RMS 
maximum deflection. The very sane results have been obtained by Dr. Ron N. 
Miles in his Ph.D. dissertation (University of Washington, 1987). 

The real challenge is the search for an analytical model that will be 
able to correlate better with the experimental data so that the strain spec- 
tral peaks are broadened 40 to 50 times instead of 2 to 3 at high SPL. The 
solutions obtained by Drs. Moyer and Miles demonstrated that time domain 
analysis with nonlinear stiffness term alone (without nonlinear damping) can 
not broaden the spectral peaks to the orders observed experimentally. Non- 
linear damping effects will have to be considered in the analysis for panels 
excited at high pressure levels. 
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(TDS - TIME DOMAIN SIMULATION) 

(MEL - METHOD OF EQUIVALENT LINEARIZATION) 
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ACTIVITY 6: Effects of Boundary Conditions and Initial Imperfections on 

Panel Response. 

To improve the agreement between predicted and measured dynamic 
response of rectangular aluninum panels, both rotational and translational 
boundary springs are considered at the panel edges. To ease the mathe- 
matical manipulations, the MACSYMA program is used. Dynamic response 
obtained using MACSYMA are compared with the improved Rayl eigh-Ri tz anal- 
ytical solutions and NASTRAN results. 

Most recently, the effects of initial imperfection and edge rotational 
restraint on large deflection random response of symmetrical ly laminated 
cross-ply rectangular plates have been studied. Figure 4 shows the RMS 
(maximum deflection/plate thickness) vs. nondimensional edge restraint K (K 
= 0 corresponds to simply supported case, K = * clamped case) at four sound 
spectrum levels from 100 to 130 dB (Re 2 x 10 N/M ) for a flat plate (W ^ 
= 0) and plates with non-dimensional initial imperfections of W ^ = 0.1, 

0.2 and 0.4. In general, initial imperfections are to increase the linear 
frequency and to reduce the RMS deflection and RMS strain (Publication P-8 
and Report R-3). 

ACTIVITY 7: Influence of Temperature on Small Deflection Vibratory Analysis 

of Panels. 

Extension of the finite element method to predict small deflection 
random response and small deflection free vibrations of thermally buckled 
rectangular plates has been completed. The formulation has been generalized 
to treat the thermal loads having nonuniform temperature distributions 
T(x,y) instead of the uniform temperature. Effects of thermal prestress, 
thermally postbuckled deformation and changes in material properties due to 
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Edge rotational restraint K 

Fig. 4. Effects of initial imperfection and edge restraint on deflection 
for a square cross-ply laminate (c = 0.01) at various SSL. 



temperature have all been considered in the formulation. Further devel- 
opment of the finite element formulation for large deflection random 
response of thermally buckled structures is continued by dames E. Locke. 
The Ph.D. dissertation committee has accepted this as his research topic on 
January 11, 1988. 
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